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,2 ' Abstract. We review recent activity searching for variations in the fundamental constants 

, of nature in quasar absorption spectra and in the laboratory. While research in this direction 

ryj ■ has been ongoing for many decades, the topic has recently been stimulated by astronomical 

O ' evidence for spatial variation in the fine-structure constant, a. This result could be confirmed 

* 5j . using different quasar data and atomic clock measurements, but there are significant challenges 

^-j' to obtain the required accuracy. We review existing measurements and discuss some of the most 

pIZh ' promising systems where any variations would be strongly enhanced. 

7— I 

, 1. Introduction 

Recently, a highly significant (~ 4<r) indication that there is a spatial gradient in values of the 
fine-structure constant, a, was reported That is, in one direction on the sky a was larger in 
the past, while in the opposite direction it seems to have been smaller. This result has massive 



• implications for the "fine-tuning" problem. It is well known that the constants of nature are 

finely tuned to allow life to exist. For example, the production of carbon from helium in stars 
(the famous triple-a reaction) is sensitive to the position of a low-energy resonance, which in 
turn is very sensitive to variations in coupling constants. If the coupling constants were slightly 
different, there would be no resonance and hence no carbon-based life. There are many other 
examples. With the detection of spatial variation of coupling constants we begin to have a 



natural explanation of fine-tuning: we simply we appeared in the region of the Universe where 
values of fundamental constants are suitable for our existence. 

We can only detect the variation of dimensionless fundamental constants. In this review we 
will discuss systems that are sensitive to variation of the fine-structure constant, a = e 2 /hc, 
the proton g- factor, g p , and the dimensionless mass ratios fi = m e /m p and X q = m q /A.QCT>- 
Here m e , m p , and m q are the electron, proton, and light-current quark masses, respectively, 
and Aqcd is the quantum chromodynamics scale, defined as the position of the Landau pole in 
the logarithm of the running strong coupling constant, a s {r) ~ 1/ln (Aqcd^/^c)- The proton 
mass m, p is proportional to Aqcd (if we neglect the few percent contribution of quark masses 
to the proton mass), therefore the relative variation of \x = m e /m p is approximately equal to 
the relative variation of X e = m e / Aqcd- In the Standard Model electron and quark masses are 



proportional to the vacuum expectation value of the Higgs field. In the same vein, g p is not a 
fundamental constant, but we can express its variation in terms of variation of light quark mass 
using the relationship Sg p /g p ~ —0.1 5X q /X q [2] (accurate calculations for nuclear g-factors can 
be found in [3]). 

The evidence for a cosmological spatial gradient in a motivates us to interpret other data 
in terms of spatial variation. However, many of the limits on variations are best interpreted as 
limits on fundamental constants other than a, since the system being examined is more sensitive 
to these other constants. A minimal hypothesis is to expect all fundamental constants to vary 
along the same direction as the a- gradient. There are some good theoretical justifications for 
this postulation. If the constants vary because they are coupled to a scalar field <3? that varies 
over space-time, for example the quintessence field <l?/c 2 or a dimensionless dilaton field, then a 
fundamental constant, X, will vary according to 

^ = k x 8$ (1) 

where kx is a dimensionless coupling coefficient. Then all constants will vary in the direction of 
the gradient of the scalar field A<£, which in turn must be directed along the axis found by the 
a- variation data in []J . 

Equation (fTJ) implies that the relative variation of different constants can be related to each 
other using equations like 

kfj, = k a (2) 

where the constants R^' can be determined from observations and compared with theories of 
spatial variation. For example, applying Grand Unification of the interactions of the Standard 
Model to the relative variation of a and a s (see, e.g., |4j|) we find 

6Xq - 35 6a 



X q a 

That is, the relative variation in X q may be much larger than the relative variation in a. The 
coefficent here is model dependent, but large values are prevalent among models in which the 
variations come from high-energy scales. Therefore, variation of fundamental constants provides 
a unique probe of the predictions of Unification Theories (see, e.g., 

In this review we discuss measurements of variation of fundamental constants in quasar 
absorption spectra and in atomic clocks. We interpret the astronomical measurements in the 
context of the observed spatial gradient in a, relating other constants to this variation using 
expressions such as ^fy. Atomic clocks may prove to be the most reliable way to corroborate 
the detection of spatial a-variation, and we discuss some of the proposed experiments where 
huge enhancements of the variation can be expected. There are many other systems which 
are not discussed here, including Big Bang Nucleosynthesis where hints of variation have been 
reported 10, 11]. We refer the reader to recent reviews (3. Il2|. 



2. Quasar absorption spectra 

We start our discussion with quasar absorption systems, since the remarkable recent detection 
of spatial variation in a from optical absorption spectra p|, will provide us with a benchmark 
for all other studies. 



2.1. Optical atomic spectra 

It is natural to analyse fine-structure intervals in the search for variation of a. Measurements 
of a-variation by comparison of cosmic and laboratory optical spectra were first performed 



by Savedoff as early as 1956 13J]. There were numerous works successfully implementing this 
"alkali-doublet" method (see review fl4|). 

In 1999, a different approach was developed: the many-multiplet method 15|, IB]- It exploits 



the fact that relativistic corrections to atomic transition frequencies can exceed the fine-structure 
interval between the excited levels by an order of magnitude (for example, an s-wave electron 
does not have the spin-orbit splitting but it has the maximal relativistic correction to energy). 
We can express the a-dependence of a transition frequency oj as 

uj = ujq + qx (3) 

where 

x = (a/a ) 2 -l«2^— (4) 

«o 

«o is the laboratory value of a and ujq is the laboratory frequency of a particular transition. 
The coefficients q can vary strongly from atom to atom and can have opposite signs in different 
transitions (for example, in s-p versus d-p transitions). Thus, any variation of a may be 
revealed by comparing different transitions in different atoms in cosmic and laboratory spectra. 
A statistical gain is also realised because many more spectral lines in different elements can be 
used. This method improves the sensitivity to any variation of a by more than an order of 
magnitude compared to the alkali-doublet method. 

Relativistic many-body calculations are used to reveal the dependence of atomic frequencies 
on a 2 (the q coefficients). We have performed accurate many-body calculations of the q 



coefficient for all transitions of astrophysical interest 15 23 . These are strong El transtions 
from the ground state in Mgl, Mgll, Till, Fel, Fell, CrII, Mil, A1II, A1III, Sill, Znll, Mnll, 
as well as many other atoms and ions which are seen in quasar absorption spectra, but have 
not yet been used in the quasar measurements because of the absence of accurate UV transition 



laboratory wavelengths. For a "shopping list" of needed measurements, see Ref. [24]. 

From the very first analyses of quasar data that utilised the new method, hints of co- 
variation were reported [25L The largest analysis, with three independent samples of data 
containing 143 absorption systems spread over redshift range 0.2 < z < 4.2, suggested that a 
was smaller in the past: 5a/ a = (—0.543 ± 0.116) x 10 -5 [22]. However, these studies all used 
spectra taken at the Keck telescope. Similar studies from another group, using our method and 
calculations but a much smaller sample of data taken at the VLT, at first showed a stringent 
null constraints (Srianand et al [28]). More careful analysis of this sample suggested that the 



errors were underestimated by a large factor [2914311] . 



The latest results, combining the Keck data and a new sample of 153 measurements from the 
VLT, shows a spatial variation in a. This gradient, which we will refer to as the "Australian 
dipole", has a declination of around —60°. This explains why the Keck data, restricted mainly 
to the northern sky since the telescope is in Hawaii at a latitude of 20° N, originally suggested 
a time-varying a that was smaller in the past. The VLT is in Chile, at latitude 25° S, giving 
the new study much more complete sky coverage. The new results are entirely consistent with 
previous ones. The "Australian dipole" of a- variation found by [l[ is 

— = (1.10 ± 0.25) x 10" 6 r cos ip (5) 
a 

where 8a/ao = (at(r) — ao)/«o is the relative variation of a at a particular place r in the 
Universe (relative to Earth at r = 0). The function rcosip describes the geometry of the spatial 
variation: ip is the angle between the direction of the measurement and the axis of the Australian 
dipole, (17.4 (0.9) h, —58 (9)°) in equatorial coordinates. The distance function is the light-travel 
distance r = ct measured in giga-lightyears. This is model dependent for large redshifts: the 



standard Acdm cosmology parametrized by WMAP5 32] is used to determine the light-travel 
time t. It is assumed here that 5a/ a^ = at zero redshift, which is supported by the data, 
however this assumption should be tested using the same absorption methods as are used at 
high redshift (e.g. by using absorbers within our own galaxy). 

The new results of a-variation are particularly striking in that if the VLT data and Keck 
data are considered separately, there is a high level of agreement between their dipole fits. In 
their paper [lj], the authors estimate the probability that the observed alignment between the 
independent Keck and VLT dipoles is due to chance at ~ 4%. Similarly, if one breaks the 
data into subsamples consisting of absorbers at z < 1.8 and those at z > 1.8 (approximately 
half the data in each subsample), then there is again agreement between the directions of the 
independent dipole axes. 

In addition to the previously mentioned many-multiplet method results of Srianand et al from 
the small VLT sample [28l ] (see also [33l]). there are some high-accuracy single-absorption system 
a-variation results using alkali-doublets (SilV) [IH and single-ion (Fell) 1 35- 3^] measurements. 
The alkali-doublet result suffers from the same problems as the small-sample many-multiplet 
result of [28l ] (namely, sharp fluctuations in chi-squared vs. 5a/ a graph which indicate failings 
in the chi-squared minimisation routine (29|). 

The Fell results give stringent constraints in two absorbing systems: at z = 1.15 towards 
HE0515-4414 [13] and at z = 1.84 towards Q1101-264 [II]. It is interesting to compare the 
results of these analyses with the variation expected in these systems if the Australian dipole 
result is correct. Using ([5]), we obtain an expected variation of 5a/a = (1.9 ± 1.5) x 1CP 6 
for the former, which was measured to be (—0.07 ± 0.84) x 10 -6 [371 ] . The absorber towards 
Q1101-264 has a larger expected variation of 5a/a = (3.8 ± 2.0) x 10~ 6 , and was measured to 
be (5.4 ±2.5) x 10~ 6 [38]. The single system results are therefore seen to be consistent with the 
Australian dipole. 

2.2. Molecular rotational quasar spectra 

Limits on variation of /i at high redshift can be obtained by comparison of different rotational- 
electronic excitations in molecular hydrogen. Taking advantage of newly available H2 
wavelengths, Ref. [3j| reported a non-zero cosmological variation of fj, in quasar spectra using 
molecular hydrogen transitions in the Ly-a forest. The authors obtained Sfx/fx = (2.4±0.6) x 10 -5 
at redshifts z ~ 2.6 — 3.0, i.e. a decrease of \i in the past 12 Gyr at the 3.5<r confidence level. 

More recent determinations, using the laboratory wavelengths and sensitivity coefficients 
presented in 3^-41] and carefully controlling systematics and accounting for known calibration 



errors in the VLT, obtain a more stringent constraint: 5fj,/fj, = (3.4 ± 2.7) x 10 [42j, |43| from 
four quasar absorption systems. 

In the context of an observed dipole in a-variation [lj], it makes sense to check whether the 



(limited) data supports interpretation as a spatial variation. In Ref. 44J the data was shown 
to support a statistically significant spatial variation of fi aligned with the Australian dipole, 
although it was noted that the paucity of data prevents any firm conclusion. 

2.3. Comparison of hydrogen hyperfine and UV transitions 

A comparison of the 21-cm hyperfine transition in atomic hydrogen with UV metal lines was 
performed for 9 quasar spectra with redshifts 0.23 < z < 2.35 45|, 4g|. The ratio of the transition 



frequencies is proportional to the parameter x = a 2 [ig p , and this was constrained to 

5x/x = (6.3 ±9.9) x 10~ 6 . (6) 



It was found that there was much more scatter in the data than would be expected from the 
statistical errors alone (46|. In principle this could mean that the model is wrong, e.g. a spatial 



Table 1. Comparison of expected variation, given by ([5]), and measured variation of fundamental 
constants in different quasar absorption systems. Rg p and are denned by equations like ([2]). 
As noted in Section [lj 1?° —O.lRg and R^ sa where i?^ is the variation of light quark 
mass X q relative to a-variation. Errors in the expected a-variation (the prefactors in the third 
column) are of the order ~ 1.5 x 10~ 6 . 
System Constant Expected variation Measurement Ref. 

(xlO" 6 ) (xl0~ 6 ) 

HI + CI cT^g~ p 1.12(2 + i£ + i2°) 6.64 ± 0.84 stat ± 6.7 sys gzT 

-5.20(2 + R« + R« p ) 7.0 ± 1.8 stat ± 6.7 sys [47] 

Hl + mol. rot. a 2 g p 0.50 (2 + R% p ) -2.0 ±4.4 [48J 

-5.47(2 + ^) -1.6 ±5.4 [48J 

HI + OH {a 2 /^- 57 g p -1.04(3.14- 1.57/2° + /%) 4.4 ± 3.6 stat ± 10 sys [49] 

OH (a 2 /fi) 1S5 g p 0.50 (3.70 - 1.85JS" + R% p ) -11.8 ±4.6 [50] 

NH 3 n -5.47 B* < 1.8 (2a) [5JJ 

1.34 Rl < 1.4 (3a) [52] 



variation should be considered rather than a time variation. However, even the best fit dipole 
model (which does not align with the Australian dipole) cannot explain the scatter [44]]. It 
is probable that velocity offsets between the hydrogen and metal absorbers cause dominating 
systematics. Two new measurements comparing HI 21-cm with UV lines in neutral carbon in 
single absorption systems have reported slightly stronger constraints in x |47J; these are shown 
in the first two lines of table [TJ 



2-4- Comparisons involving hyperfine and molecular rotational transitions 

Individual measurements of fundamental constant variation in quasar absorption spectra can 
be compared with the observed spatial a-variation [3]. That is, we can regard the a-variation, 
according to ([5]) with r and ip calculated for the object, as the "expected" variation and compare 
other measurements with it (this is presented in table [1]). 

The frequency of the hydrogenic hyperfine line is proportional to a 2 fig p ; molecular rotational 
frequencies are prop ortional to fx. Comparison of the two placed limits on variation of the 
parameter a 2 g p [53J in two quasar absorption spectra. A similar analysis was repeated using 
more accurate data for the same objects [48], at z = 0.247 and at z = 0.6847, resulting in the 
limits shown in lines 3 and 4 of table [TJ respectively. The object at z = 0.6847 is associated with 
the gravitational lens toward quasar B0218±357 and corresponds to lookback time ~ 6.2 Gyr. 

Comparison of OH 18-cm with H I 21-cm, and different conjugate-satellite OH 18-cm lines, can 
be used to measure the combinations of fundamental constants (a 2 ///) 1 ' 5 g p and (a 2 / /i) 1 ' 85 g p , 
respectively. Resulting measurements are shown in lines 5 and 6 of table [TJ 



2. 5. Enhancement of variation of // in the inversion spectrum of ammonia 

In 2004, van Veldhoven et al suggested usin g a decelerated molecular beam of ND3 to search for 

the variation of \i in laboratory experiments [54]. The ammonia molecule has a pyramidal shape 

and the inversion frequency depends on the exponentially small tunneling of three hydrogen 

(or deuterium) nuclei through the potential barrier. Because of that, it is very sensitive to any 

changes of the parameters of the system, and particularly to the reduced mass for this vibrational 

mode. 

High precision data on the redshifts of NH3 inversion lines exist for the previously mentioned 
object B0218±357 at z = 0.6847 55J]. Comparing them with the redshifts of rotational lines of 
CO, HCO + , and HCN molecules from Ref. [56J allows strong limits on variation of fi, and hence 



X e , to be obtained 57J. The most accurate measurements, utilising new rotational spectra, give 
very stringent limits on /x-variation in this quasar absorption system and the object PKS1830- 
211 at z = 0.8858 (last two lines of tabled]). 



3. Clocks 

3.1. Optical atomic clocks 

Atomic clocks can be used to measure time-variation of fundamental constants in the Earth 
frame. Different optical atomic clocks utilize transitions that have positive, negative or small 
contributions of the relativistic corrections to frequencies, so comparison of these clocks can be 
used to measure a-variation. The same methods of relativistic many-body calculations used 
in the quasar absorption studies can be used to calculate the dependence on a of different 



clocks [la . 1 1 71 . l58H6ll|. A summary of results is presented in Ref. [62]]. The relativistic effects 
are proportional to (Za) 2 , therefore the q coefficients for optical clock transitions may be 
substantially larger than in cosmic transitions since the clock transitions are often in heavy 
atoms (Hgll, Ybll, Yblll, etc.) while cosmic spectra contain mostly light atoms (Z < 33). 

The temporal variation of a measured in the laboratory can be compared with the variation 
expected from the observed spatial gradient in a [l[ because the solar system moves along the 



axis of the dipole. The expected variation from the dipole was calculated to be [63J] 

a/a = 1.35 x 1(T 18 cos ip yr" 1 (7) 

where ip is the angle between the motion of the Sun and the dipole. The best fit from [l[ gives 
costp ~ 0.07 but this has an uncertainty of ~ 0.15. The a- variation is modulated by the annual 
motion of the Earth around the Sun, 5a/ a = 1.4 x 10~ 20 cosu;i [63] where u refers to the angular 
frequency of the yearly orbit. 

The current best constraint on time-variation of a was achieved by precisely measuring the 



frequency ratio of Hg II and Al II clocks several times over the course of a year 64| . Using our 
calculations, the rate of change of a is measured at a/ a = (—1.6 ± 2.3) x 10~ 17 yr _1 . This limit 
will need to be improved by around two orders-of-magnitude in order to confirm or contradict 
the observed spatial gradient in a. On the other hand, the quasar observations do not exclude 
time- variation of a below the rate of ~ 10~ 16 yr" 1 , and here the current laboratory limits are 
already competitive. 



3.2. Enhanced effect of a-variation in Dy atom 

The sensitivity required by © may be obtained by finding systems where a-variation is strongly 
enhanced. Transitions between two almost degenerate levels in Dy atom can give a very high 



15, 17, mi 



relative enhancement because these levels move in opposite directions if a varies 
The relative variation may be presented as 5u/lo = K5a/a where the coefficient K exceeds 10 s 
(q = 30 000 cm -1 , uj ~ 10~ 4 cm -1 ). The values of K = 2q/co are different for different hyperfine 
components and isotopes because u changes. An experiment is currently underway to place 
limits on a variation using this transition [6£ 



661 ] . however one of the levels has quite a large 



linewidth and this limits the accuracy. The current limit is a/a = (—2.7 ± 2.6) x 
Several other enhanced effects of a variation in atoms have been calculated [67J, K> 



10 



-15 



-1 



3.3. Enhanced effects of a-variation in highly charged ions 

Sensitivity to a-variation increases with ion charge as {Z{ + l) 2 The most sensitive 

atomic systems will maximize the contributions from three factors: high nuclear charge Z, 
high ionization degree, and significant differences in the configuration composition of the states 
involved. Unfortunately, the interval between different energy levels in an ion also increases 
as ~ (Zi + l) 2 , which can quickly take the transition frequency out of the range of lasers 



as Zi increases. The phenomena of Coulomb degeneracy and configuration crossing can be 
used to combat this tendency. In a neutral atom, an electron orbital with a larger angular 
momentum is significantly higher than one with smaller angular momentum but with the same 
principal quantum number n. On the other hand, in the hydrogen-like limit orbitals with 
different angular momentum but the same principal quantum number are nearly degenerate. 
Therefore, somewhere in between there can be a crossing point where two levels with different 
angular momentum and principal quantum number can come close together: in such cases the 
excitation energy may be within laser range. 

In Ref. [6^] we showed, using the Ag isoelectronic sequence as an example, why high q- 
values can occur in highly charged ions, and how the tendency of such systems towards large 
transition frequencies could be overcome. A two- valence-electron ion, Sm + , was identified, 
which has optical transitions that are the most sensitive to potential variation of a ever found. 
While atomic spectroscopy in electron beam ion traps is currently not competitive with optical 



frequency standards (see, e.g., [70j, |71[ and review [72j|) the technology continues to improve, 
and with the enhancements in sensitivity, highly-charged ions may prove to be a good system 
for detecting variation of a. 

3.4- Enhanced effect of variation in UV transition of 229 Th nucleus 

The 229 Th nucleus has the lowest known excited state, lying just 7.6 ± 0.5 eV above the ground 
state [73]. The position of this level was determined from the energy differences of many high- 
energy 7-transitions to the ground and first-excited states. The subtraction produces the large 
uncertainty in the position of the 7.6 eV excited state. The width of this level is estimated to 



be about 10 -4 Hz 74], which explains why it is so hard to find the direct radiation in this very 
weak transition. Nevertheless, the search for the direct radiation continues. 

Because the 229 Th transition is very narrow and can be investigated with laser spectroscopy, 
it is a possible reference for an optical clock of very high accuracy [75J. The near degeneracy 
of these isomers is a result of cancellation between very large energy contributions (order of 
MeV). Since these contributions would have different dependences on fundamental constants, 
this transition would be a very sensitive probe of possible variation of fundamental constants [76T ]. 
A rough estimate for the relative variation of the 229 Th transition frequency is 

5u tr / 5a 5X a \ , , 

— ~ 10 5 0.1— + — 1 . 8 
oj \ a X q J 

Therefore, the experiment would have the potential of improving the sensitivity to temporal 
variation of the fundamental constants by many orders of magnitude. 

More accurate nuclear calculations give different values for the sensitivity of this transition 



to a. Refs. [771 . [78J claim that both isomers have identical deformations and therefore there is 
no enhancement of a-variation. Other calculations giv e enhancement factors in the range 10 2 
- 10 5 , depending on particulars of the model used [7914821] . To resolve this, we have proposed 



a method of extracting sensitivity to a-variation using direct laboratory measurements of the 
change in nuclear mean-square charge radius between the isomers (isomeric shift) [83]. 
From (|8|), we obtain the following energy shift in the 7.6 eV 229 Th transition: 

ft X 

iww-^MeV. (9) 

X q 

This corresponds to the frequency shift <5zy ~ 3 • 10 20 5X q /X q Hz. The width of this transition 
is 10 Hz so one may hope to get the sensitivity to the variation of X q about 10 -24 per year. 
This is 10 10 times better than the current atomic clock limit on the variation of X q . 



Note that there are other narrow low-energy levels in nuclei, for example the 76 eV level 
in 235 U with lifetime 26.6 minutes is the second-lowest known. One may expect a similar 
enhancement there. Unfortunately, this level cannot be reached with usual lasers. In principle, 
it may be investigated using a free-electron laser or synchrotron radiation. However, the accuracy 
of the frequency measurements is much lower in this case. 
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